Generalized persistence exponents: an exactly soluble model 
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It was recently realized that the persistence exponent appearing in the dynamics of nonequilib- 
rium systems is a special member of a continuously varying family of exponents, describing gener- 
alized persistence properties. We propose and solve a simplified model of coarsening, where time 
intervals between spin flips are independent, and distributed according to a Levy law. Both the limit 
distribution of the mean magnetization and the generalized persistence exponents are obtained ex- 
actly. 
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The surprise caused by the discovery of new nontriv- 
ial exponents in the dynamics of simple nonequilibrium 
systems motivated a long series of works, mainly 
devoted to the search of simple models or experimen- 
tal situations, where the so called persistence exponents 
could be computed or measured [p|"p^|. More recently, 
the idea of persistent large deviations |T^] led to the in- 
troduction of families of new nontrivial persistence expo- 
nents in, e.g., the one dimensional Glauber- Ising chain at 
zero temperature, or the simple diffusion equation. The 
probability of persistent deviations, defined as the prob- 
ability that the mean magnetization Alt = t~^ dua(u) 
of the spin (or of the sign of the diffusing field) at a given 
site was, for all previous times, greater than some level x, 
where — 1 < a; < 1, was observed to decay algebraically 
at large times, with an exponent 0{x) continuously vary- 
ing with X. When x — 1, this probability is the usual 
persistence probability, since imposing Mt = 1 requires 
that the spin never flipped. Hence 6{1) = 6, the usual 
persistence exponent. Furthermore, the distribution of 
Mt does not peak around zero for t ^ oo, but tends to a 
nontrivial limit distribution on [—1,1], singular at both 
ends as (1 =f x)^^-^. For coarsening systems, computing 
the exact value of 9 turns out to be a hard problem, so 
one does not expect the computation of 0{x), or even of 
the distribution of the mean magnetization, to be easily 
reachable. The origin of the difficulty is that spins at 
different sites are strongly correlated. 

The aim of this letter is to present an extremely sim- 
plified version of the coarsening models mentioned above, 
which allows for exact analytical expressions both of the 
limit distribution of the mean magnetization M , and of 
the generalized persistence exponents 6{x). Despite its 
simplicity, this model retains the essential features of the 
coarsening process, in particular its non stationary prop- 
erties, as will be discussed below. 

In this model, which describes the dynamics of a sin- 
gle spin, we assume that the time intervals between spin 
flips are independent. It is indeed intuitively clear that. 



because of the ever growing size of domains in coarsen- 
ing systems (or of the diffusion length in the diffusion 
equation), a spin at a given site can remain in the same 
direction for a very long time before a domain wall crosses 
this particular point and flips the spin in the reversed di- 
rection. By definition of the persistence exponent 9, the 
time r before a spin is flipped is very broadly distributed, 
with a power law tail decaying as r"^~^ for large r. The 
simplest approximation is therefore to neglect the corre- 
lations between the different time intervals between flips, 
all assumed to be distributed with the same density p{t), 
decaying as r~^~^. 

For simplicity, the distribution of time intervals p{t) 
is chosen to be a positive stable Levy distribution of in- 
dex < 6* < 1 denoted by L^(t). (The case > 1 
will be discussed below.) Its Laplace transform reads 
Lg{s) — exp(— &s^), where b is the scale factor of the dis- 
tribution, i.e. the typical values of r are of order 
p2| . As is well known, ^^(t) decays asymptotically as 
r^^^ Q. We will always suppose that a{t = 0) ~ 1. 
On the time axis, the process thus defined is a renewal 
process. 

We have investigated the statistics of the process, both 
after n sign changes, or at time t, with very similar re- 
sults in the asymptotic regime. After n sign changes, the 
time elapsed and the magnetization of the spin read 
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while, at time t, they are given by 

t = tN,+\, St = tMt^ Sn, + {-)^'\. (3) 

In the first case, n is given and tn is a random variable, 
while in the second one, t is given and Nt is the random 
variable equal to the largest n for which t„ < t. Finally 
A is the length of time measured backwards from t to the 
last event. The corresponding distributions are defined 
as 
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P(n, x)=P (M„ = Sn/tn > x) , 
P{t,x) ^P{Mt ^ St/t>x). 



(4) 
(5) 



For distributions which are peaked around their means at 
large times, these quantities are referred to as the proba- 
bihties of large deviations and are exponentially decreas- 
ing with n or i respectively. In the present case, where 
p{t) = Lg^r) is a positive Levy distribution, we find the 
limit distribution 



P{x) — lim P{n,x) — lim P{t,x), 

n — >oc t — *oo 
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where uj = (1 — a;)/(l + x) and r = 1 (see below). 



(6) 
(7) 



10" 



R(t,x) 


R(n,x) 



















10" 



10 ' 



- 10" 



- 10"' 



FIG. 1. Plot oiR{t,x) (left) and R{n,x) (right) for 6 = 1/2 
and various values of x, in log-log coordinates. The power-law 
behavior of both quantities for large times is clearly seen. 

Let us sketch the proof of (0) for P{n,x), leaving the 
calculation of P{t, x) to a longer publication jl^. We in- 
troduce and , which are the lengths of time spent 
by the spin, respectively in the positive or negative di- 
rection, such that tn = T+ -I- T~ and Sn ~ — , 
with 1"+ = Ti + T3 + • • • + T2fe+i, if n = 2fc + 1, and T+ = 

n+TsH VT2k-i, if n = 2fc, andr^ = T2+Ti^ VT2k, 

in both cases. Then P{Sn/tn > x) ^ P{T- /T+ < uj) 
with UJ — {1 ~ x)/{l + x). Since and are sums 
of stable Levy random variables Lg, they are themselves 
stable Levy random variables Lg , where, using the ad- 
dition rule of the scale parameters, b~ — kb, and b'^ = kb 
(if n = 2fc), or b+ = {k + 1)6 (if n ^ 2k + 1). The deter- 
mination of P{n, x) therefore amounts to computing the 
distribution of the ratio of two Levy laws with parame- 
ters b~ and fe"*". Denoting by H the Heaviside function, 
and using its Laplace representation along the Bromwich 
contour, one finds 



p(T-/r+ > u) 

= / ^exp[-6+(sc^)Vxp[-6"( 



This integral leads to (|^) with r = b~ /b~^ . In the limit 
n oo, r 1. This derivation also shows that whenever 
n is even, P{n, x) ~ P{x). 

The limit density f{x) = —P'ix) of the mean magne- 
tization reads 
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It is even, and diverges when x ±1 as (1 =F x)^^^. 
For 6 < 9c = 0.5946 . . ., where 6c is the solution of 
6c — cos {tt9c/2), X — corresponds to a minimum of 
/(x), while for larger 6, it corresponds to a local max- 
imum. This can be interpreted as a precursory sign of 
the fact that f{x) tends to 5{x) for 6 > 1. (It also shows 
that f{x) cannot be approximated by a beta distribu- 
tion when 6 is too large. In this respect, compare to the 
discussion in |Q.) 
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FIG. 2. Plot of the exponents (^(x) and B(x)ie for 
d = 0.3, 0.5, 0.7, showing that the relation Q(x) = 6(j}{x) holds. 
The lines corresponds to the exact result ^(a::) = 1 — P{x). 



We now consider the probability of persistent large de- 
viations, defined as the probability that the mean mag- 
netization M was, for all previous times, greater than 
some level x. More precisely one defines the quanti- 
ties R{n,x) — P(M„' > x,Vn' < n) and similarly 
R{t,x) = P{Mt' > x,yt' < t). Numerical computa- 
tions show that both quantities decay algebraically in 
the asymptotic regime (see Fig. 1), respectively as 

R{n, x) - n-'^(^) {n > 1), R{t, x) - t^^f^' {t > 1), 

where the two families of exponents are related by 6{x) = 
6<j){x) (see Fig. 2). This relation is indeed expected since 
for a given n, tn scales as n^^^ . Note that by definition 
of the model, 6{1) = 6. We also observe with very good 
accuracy (see Fig. 2) the relation 



(x) = 1-P{x) = j duf{u), 



(9) 
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which we now estabhsh exactly. 

For this, we note that R{n, x) is the joint probabil- 
ity that Sn' > xtn' for all 1 < n' < n. Since clearly 
R{2k,x) — R{2k + l,x), we assume that n is even, and 
write 

R{n - 2k, x) = P(ei> 0, 6 + 6 > 0, . . . , 

ei+6 + ---+a->o), (10) 

where = (1 ~ x)T2i~i — (1 -I- x)T2i- Since the 
are positive Levy variables of index 9, the are also 
Levy variables of index 9, with an asymmetry parame- 
ter (3 = {bj^ — l)/(tf^ + l)j which measures the relative 
weight of the negative and positive tails . The solu- 
tion to (^o|) for general stable Levy variables is known 
io|,f5|. It reads 



R{n = 2k, x) 



r(fc + 1 - g) 

T{k + l)T{l-qy 



(11) 



where 1—q is the probability that ^ is positive. This prob- 
ability is precisely the quantity P{n — 2,x) introduced 
above, itself equal to P{x). Hence q = 1 — P{x). Finally, 
the large k behavior of the r.h.s. of ( pi] ) is cx i.e., 
4>{x) = q, which completes the proof of (^. Note that 
the plot of ( |Tl| ) is indistinguishable from that obtained 
numerically for R{n,x). Eqs. (Q), (||), and (^l]) are the 
main results of this work. 

The rest of this paper is devoted to a discussion of 
further properties of the model, and to a presentation of 
some possible generalizations. 

First, the stochastic process presented above, where 
time intervals between spin flips are independent and 
distributed according to a Levy distribution, exhibits 
nontrivial temporal prop erties, both from mathematical 
|^,|2|, and physical [^Bp^s) points of view. For example, 
although p(r) is fixed in time, the probability distribution 
of the length of time A from some time origin (or wait- 
ing time) tw to the next flip is non stationary for 9 < I, 
i.e. it depends both on and A, while it is asymptoti- 
cally independent of for 9 > 1. As a consequence, the 
probability that a given spin did not flip between times 
iw and + t is a function of ty^/t if 9 < 1, while it is 



16[. Thus for 9 < 1, this 
nature of the persistence 



independent of if 9 > 1 
model captures the aging [09 
phenomenon. This property is deeply related to the fact 
that the largest Ti in the sum t„ = X]"=i contributes to 
a finite fraction of tn for 9 < 1 even in the limit n — s- cx), 
while this fraction is asymptotically zero for 9 > 1 ||l^,|l8| . 
Correspondingly, this also ensures that the distribution 
of the mean magnetization does not peak around x = 0, 
as was shown above. 

Despite its simplicity, the model discussed here thus 
shares many features of more complex coarsening pro- 
cesses. As shown above, it leads to nontrivial predic- 
tions for the quantities P{x) and 9{x). Also, the be- 
havior of R{t, x) observed in Fig. 1 strongly resembles 



that found in for the Glauber-Ising chain or the dif- 
fusion equation. These predictions can be seen as ap- 
proximations for these more general models. In Fig. 3, 
we compare, for the Glauber model at zero temperature, 
the function 9q{x)/9g{^), as determined numerically in 
pl| , both with 1 — P{x), where P{x) is given by (0) with 
6* = 3/8 1^, and with 1 — Pg{x), the distribution of mag- 
netization measured numerically in Jill . Although there 
is qualitative agreement between the three curves, the 
above relations are clearly only approximate. Further- 
more there remains to understand the qualitative differ- 
ence in behavior between the persistence exponents 9{x) 
for diffusion and for the Glauber-Ising chain. It would 
be therefore interesting to generalize the present model 
to include some correlations between the time intervals 
Ti. The independent interval model presented here is ac- 
tually, in many respects, similar to the random energy 
model (rem) for spin-glasses |17 1^. For example, the 
r.h.s. of (^l]) is identical to the expression for the partici- 
pation ratio Yfe+i in the rem, with a reduced temperature 
equal to q 0Jl|,|l|l. 
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FIG. 3. Comparison between the function SG(a^)/^G(l) for 
the Glauber-Ising model in one dimension, as determined 
numerically in and two predictions inspired from the 

present model: 1 — P{x), where P{x) is given by (^) with 
e = 3/8, and 1 - Pcix) (see [0). 

The only feature of p{t) relevant for the results given 
here is its asymptotic power-law behavior at large values 
of n or t, and not its detailed behavior for small r. The 
reason is that the sum of a large number of power-law dis- 
tributed variables converges (for 9 < 1) towards the Levy 
distributions considered here. In our simulations we used 
two different distributions both decaying as r^^~^ for 
large t, obtaining the same results for sufficiently large 
times. Note that the relations (l){x) ~ 9{x)/9 = 1 — P{x) 
are best obeyed numerically for 9 = 1/2, because in this 
case it is easy to generate the corresponding stable dis- 
tribution. 

We also studied the case 6* > 1, where p(r) has a 
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finite first moment. In this case, it is easy to check 
that tn grows Unearly with n, while Sn grows as n^^^ 
for 1 < 6* < 2 and as ^/n for 9 > 2 Hence, the 

quantity A/„ tends to zero for large n, and f{x) col- 
lapses to a 6 function. However, the persistence expo- 
nents 9{x) remain well defined, and are found to be equal 
to 0{x > 0) = 6*, e{x = 0) = 1/2 and e{x < 0) = 0. This 
shows that the relation between 9{x) and P{x) actually 
still holds in this degenerate case, except for a: = where 
the value of P{x) is ill defined. However, the nature of 
the persistence phenomenon in this model is quite differ- 
ent when 6 > 1, where it becomes stationary (see above). 
It would be interesting to see if this is also true of more 
general models where 9 > 1, such as the diffusion equa- 
tion in high dimensions P, |lT| , ^ . 

We have generalized slightly the problem, by choosing 
the time intervals during which the spin cr is respectively 
equal to 1 or to —1 with a different scale factor. The 
distribution f{x) becomes asymmetrical. One can how- 
ever check, both numerically and analytically, that the 
relations 9{x)/9 — (t>{x) = q = 1 — P{x) still hold in this 
case. We have also relaxed the deterministic alternation 
of signs, and considered Sn = X^iLi where the a.; 
are independent identically distributed random variables, 
with (a) = and {a?) finite. The above results (|^) and 
for the limit distribution and the exponents remain 
unchanged. 

When 9 = 1/2, p{t) = L'^^^{t) is precisely the distribu- 
tion of the time intervals between two returns to the ori- 
gin of the binomial random walk with equal steps ±1, in 
the regime of long times. In this sense the binomial ran- 
dom walk is 'primitive' with respect to the present 'walk' 
with time-space coordinates t, St (or Sn), and instan- 
taneous velocity equal to cr(i). For the latter, y = /t, 
is the fraction of time spent by the walk stepping in the 
positive direction, or for the primitive walk, the fraction 
of time spent on the positive half axis. Its distribution 
is well known, and given at large times by the arc sine 
density \/^:^Jy{l — y), which is precisely the result (||), 
with 9 = 1/2, and x = St/t = 2r+/i- 1. In this respect, 
(^ appears as a generalization of the arc sine law. A 
striking consequence of the present work is the existence 
of the families of exponents 9{x) and 4>{x), since, when 
9 = 1/2, these exponents describe properties of the sim- 
ple random walk. We note that the present work provides 
a clue to the determination of the temporal behavior of 
the hierarchy of quantities introduced in [Tl|] . 

Finally, it is tempting to conjecture that for a generic 
stochastic process St such that there exist an a > for 
which St /i" admits a nontrivial limit distribution at large 
times, then 

p(f^>2;,W'<t)(xt-''(-) (12) 

Such a conjecture, verified in our model, is also corrob- 
orated by the works of pl|,p2| on the random walk (for 



which a = 1/2). In this case x is not restricted to a finite 
interval, and 9{x) cannot be simply expressed in terms 
of the Gaussian distribution of St/Vi- 

We wish to thank J.M. Luck, S. Redner, and C. Sire 
for interesting discussions. 
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